rho = [
    0.1
0.166666667
0.2
0.2
0.25
0.285714286
0.32
0.366666667
0.4
0.45
0.7
0.909090909
    ];

L_teorico = [
0.111111111
0.2
0.25
0.25
0.333333333
0.4
0.470588235
0.578947368
0.666666667
0.818181818
2.333333333
10
    

];

L_medido = [
    2.383
2.558
2.698
2.698
2.951
3.061
3.324
3.469
3.677
3.959
6.258
7.424

];

hold on
plot(rho, L_teorico, ':k');
plot(rho, L_medido, '-k');
